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Temperature maps of the Cosmic Microwave Background (CMB) radiation, as those obtained by 
the Wilkinson Microwave Anisotropy Probe (WMAP), provide one of the most precise data sets to 
test fundamental hypotheses of modern cosmology. One of these issues is related to the statistical 
properties of the CMB temperature fluctuations, which would have been produced by Gaussian 
random density fluctuations when matter and radiation were in thermal equilibrium in the early 
Universe. We analysed here the WMAP data and found that the distribution of the CMB temper¬ 
ature fluctuations P CMB (AT) can be quite well fitted by the anomalous temperature distribution 
emerging within nonextensive statistical mechanics. This theory is based on the nonextensive en¬ 
tropy S q = fc{l — f dx[P q (x)] q }/(q — 1), with the Boltzmann-Gibbs expression as the limit case 
q — > 1. For the frequencies investigated ( v = 40.7, 60.8, and 93.5 GHz), we found that P CMB (AT) 
is well described by P q (AT) oc 1/[1 + ( q — l)B(y)AT 2 ] 1 ^ ,_1 ', with q = 1.055±0.002, which exclude, 
at the 99% confidence level, exact Gaussian temperature distributions p Gauss (AT) oc e - B (") AT , 
corresponding to the q —► 1 limit, to properly represent the CMB temperature fluctuations measured 
by WMAP. 
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Cosmic Microwave Background (CMB) radiation is formed 
by the leftover photons from the early Universe, when mat¬ 
ter and radiation were coupled in thermal equilibrium. With 
the expansion of the Universe, these photons decoupled and 
spread out freely throughout space, basically conserving their 
primordial features. In the early 90’s, the Far Infrared Ab¬ 
solute Spectrophotometer, on board the Cosmic Background 
Explorer (COBE) satellite, proved that this radiation was in 
thermal equilibrium by measuring its Planckian spectrum, 
with the temperature To = 2.725 ± 0.002 K (lj. Although 
the accuracy of these data (within limits as tight as 0.03% 
in the frequency range 60-600 GHz) left no doubt about the 
past thermal equlibrium state of the CMB, it is still possible 
that such Planck law derived within Boltzmann-Gibbs statis¬ 
tics does not exactly describe this radiation and may instead 
obey a generalized expression. A plausible distribution could 
be that obtained within nonextensive statistical mechanics 
(for a review on this subject see, e.g. Q), with values of the 
Qa-parameter ranging in the interval \q a — 1| < 5 x 10~ 4 (a 
stands for average and refers to the quasi-Planckian distribu¬ 
tion corresponding to To) as shown in ,3:]. 

Observations with another instrument on board COBE, 
the Differential Microwave Radiometer (COBE-DMR) 0, de¬ 
tected for the first time that the CMB contains tiny variations 
around To, termed CMB temperature fluctuations AT, at the 
level of one part in 10 5 on large angular scales (~ 7°). Since 
the standard inflationary cosmology predicts that these tem¬ 
perature fluctuations should be isotropic and Gaussian ran¬ 
dom, the COBE-DMR data motivated a number of analyses, 
although not so accurate due to the large angular resolution 
of the data, to test the statistical properties of the CMB (see, 
e.g. 0). 

Recently, highly precise and excellent angular resolu¬ 


tion data from the Wilkinson Microwave Anisotropy Probe 
(WMAP) fe] confirmed the existence of the CMB temperature 
fluctuations. The WMAP satellite observes the microwave 
sky in five frequency bands, K, Ka, Q, V, and W, centred 
on the frequencies 22.8, 33.0, 40.7, 60.8, and 93.5 GHz, re¬ 
spectively. The corresponding CMB maps released by the 
WMAP team are pixelized in the HEALPix scheme 0] with 
a resolution parameter N s id e = 512, which means that the 
celestial sphere is covered with 3,145, 728 equal-area pixels, 
with a pixel size of ~ 7 arcmin. 

These highly accurate CMB data have renewed concerns 
over the CMB statistical properties, and considerable analy¬ 
ses of the Gaussian hypothesis have been done f(|. Clearly, 
the study of such hypothesis must take into account the 
possibility that deviations from Gaussianity may have non- 
cosmological origins such as unsubtracted foreground con¬ 
tamination, instrumental noise, and/or systematic effects §Jj. 
But they may also have cosmological origin, as being, for in¬ 
stance, the effect of cosmic strings on CMB M- In a variety 
of analyses using different mathematical tools, and includ¬ 
ing foreground cleaning processes aimed to eliminate possible 
non-Gaussian contaminations, many evidences regarding de¬ 
viations from Gaussianity in the WMAP CMB data have been 
recently reported na (see also HI and references therein). 

In what follows we shall perform the statistical analysis 
of the Q, V and W maps, after the application of the KpO 
mask to eliminate known foregrounds, in order to determine 
how much their distributions of temperature fluctuations de¬ 
viate, if they do, from the Gaussian distribution. Then we dis¬ 
cuss the possible account of such distributions according to 
the non-Gaussian temperature distribution emerging within 
nonextensive statistical mechanics, because gravitation is a 
long-range interaction and this kind of phenomenum seems to 
be conveniently studied in the framework of this theory 0 . 
Definitely, the statistical significance of our results shall be 
supported by the analysis of substantial Monte Carlo CMB 


2 


maps. 

Now, we briefly introduce the basics of nonextensive statis¬ 
tical mechanics (see III for various applications). The prob¬ 
ability distribution P q (x), as a function of the variable x, re¬ 
sults from the optimization of the g-entropy defined by [if 


S q = k{l- J dx [P q {x)] q }/{q - 1), 

(1) 

with the constraints 


J dx Pq (x) = 1 , 

(2) 

f dxx [P q (x)] q 

Jdx[P q {x)]i 

(3) 

f dxx 2 [Pq(x)] q 

J dx[Pq(x)]i 

(4) 


From this, one straightforwardly obtains 


P*(x) = 


- Xlq{x-C 1 )-X 2 q(x 2 -C 2 ) 


fdx 


-\l q (x-ci)-\2 q (x 2 -c 2 ) ’ 


(5) 


where Xi q and \ 2 q are related to the Lagrange multipliers, 
and where 


e q = l 1 + (! - <l)z] 1,{1 9 \ for [1 + (1 — q)z] > 0 , (6) 

while e q = 0 otherwise. Note that this distribution function 
P q is the solution of a non-linear Fokker-Planck equation 0- 
One observes that eq. © can be rewritten as follows 


are removed from the beginning with the application of the 
KpO mask 0 . Thus, one is left with CMB signal plus the 
Gaussian signal from the instrument noise (actually the signal 
noise is Gaussian per observation ii). 

Thus, one has to consider these effects by defining |l0| the 
variance of the total signal Tpj xe j as 

2 

2 _ 2 a 0 . 2 ,\ 
a » = °pixel - — + °CMB , (11) 

where rii is the number of observations for the ith pixel, ag is 
the noise variance per observation, and cr^ MB is the variance of 
the CMB temperature fluctuations. The mean contribution 
of the instrumental noise can be estimated by considering the 
effect of the different number of observations for each pixel 
(see Ref. [l|j, pag. 16). Thus, given a CMB map, this can be 
done with the effective variance due to noise 
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°noise 




^o 2 Eti(l /m) 
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( 12 ) 


where ag is the variance per observation characteristic of the 
instrument (radiometer), n; is the effective number of ob¬ 
servations for the ith pixel, and N is the total number of 
pixels considered in the analysis of the map. For the Q, V 
and W maps we obtained E(l/ n ») = 7, 423.90, 5, 465.98, and 
1,816.96, respectively; moreover the number of pixels anal¬ 
ysed N is equal for the three maps Nq = Ny = Nw = 
2,414, 705. In other words, the effective noise variances 
°noise M an< ^ ^he variances crjj leads to the CMB variance 


P q {x) 


-B q (x-Xo ) 2 

Cq 


f dxe 


— Bq(x — x o) 

q 
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-B q (i-i 0 ) 2 

c-q , 


(7) 


where the values of B q and xo are related to Ai 9 and X 2 q ; 
A q is the normalization constant obtained in such a way that 
eq. © is satisfied. 

We shall now use this nonextensive distribution to account 
for the CMB temperature fluctuations distribution of WMAP 
data. In this case we have x = T, x o = To, and B q = B(y) 
(strictly speaking B q = B q (y)\ however, for the cases consid¬ 
ered here we have just one value for the q-parameter). Thus, 


P q (AT) = A q e q BMAT2 , (8) 

where we wrote P q (AT) instead of P q (T) to be clear that 
our analysis is dealing with the statistics of the temperature 
fluctuations. In the limit q —> 1, we recover the Gaussian 
distribution 

Pq ^ p Gauaa = A e -B(„)AT 2 ^ (Q) 


where A = l/(er„V27r), B(y) = 1/(2crj), and is the vari¬ 
ance of the Gaussian distribution. For a comparison of the dis¬ 
tributions P q with p Gauaa i i n our analyses we assume A q = A. 

As recently pointed out by Jeong & Smoot 0, the signal 
measured at any pixel in the microwave sky is made of several 
components 


Tpixel 7f ore g rounc i s + T no j se + T CMB , (10) 

corresponding to foreground signals, the noise from the in¬ 
struments, and the CMB temperature fluctuations, respec¬ 
tively. Foreground contributions, expected to be small away 
from the Galactic plane and with point sources punched out, 


CCMB = - cr noise( J/ ) ■ ( 13 ) 

Although both al and o'noiseM depend on the map under 
analyses, their difference is independent of the map, in other 
words, if our treatment of instrumental noise is correct, the 
CMB variance o"cmb should be the same for the three CMB 
maps under investigation (Q, V, and W). 

As previously reported (see section 2.2 and Fig. 2 in 
Ref. [loj) the WMAP CMB temperature distribution does 
not fully obey a Gaussian temperature distribution 


S (AT) = 


y/2n 


0 -(l/2crJ)AT 2 


(14) 


In fact, the deviation from a Gaussian distribution can be 
apreciated in Figs, la, lc, and le, for the Q, V, and W maps, 
respectively. In our analyses, the best-fit Gaussian to the data 
was obtained according to the x 2 /degree of freedom (dof) es¬ 
timator test. Thus, the x 2 /200 values for the Gaussian fits 
are 0.116, 0.275, and 0.167, for the Q, V and W maps, respec¬ 
tively, and the x 2 /200 values for the nonextensive tempera¬ 
ture distribution P q are 0.00155, 0.00178, and 0.00216, for the 
Q, V and W maps, respectively. Analyses performed using 
400 dof instead of 200 dof result in similar estimative values. 

In Fig. 1 the blue lines correspond to the best-fit Gaus¬ 
sian with variances a q = 140.58/rK, av = 163.52/rK, and 
<7w = 190.35/rK, while the red lines correspond to the best 
nonextensive distribution fit with the same variances and 
q = 1.056. Once the variances have been determined 
through the \ 2 best-fit Gaussian temperature distribution, 
for each of the CMB maps, then we use the effective noise 
variance cr^ 0 j se given in eq. o again for each of the CMB 
maps, to calculate the CMB variance. Our results are oq MB = 








3 


(68.77) 2 , (69.34) 2 , (68.81) 2 /rK 2 , for the Q, V, and W maps, re¬ 
spectively, in excellent agreement with what is expected. This 
result validates the effective noise variance as representing the 
mean contribution of the instrumental noise. 

In Figs. la, lc, and le we plotted the temperature 
distributions of the Q, V, and W maps, in the form 
logio(P CMB data ),log 10 (P Gauss ), and log 10 (P 9 ) versus AT us¬ 
ing the x 2 best-fit Gaussian temperature distribution p Gauaa 
and also the best-fit nonextensive distributions P q with 
q = 1.056. To enhance the non-Gaussian behavior of the 
WMAP data, in Figs. lb, Id, and If, we plotted in¬ 
stead log 10 (P CMB data ),log 10 (P GauBB ), and log 10 (Pj) versus 
(AT/ov) 2 , since linearity (blue curves) corresponds to the 
best-fit Gaussian distribution. Our results evidence that the 
distibution of the CMB temperature fluctuations does not 
obey a Gaussian distribution. 

In order to strengthen our analysis, we removed fits to 
the WMAP foreground templates for the eight individual ra¬ 
diometer maps (corresponding to radiometers Ql, Q2, VI, 
V2, Wl, W2, W3, and W4), after applying the KpO mask to 
avoid high latitude foregrounds. Our results are plotted in 
Fig. 2. As observed, these results were essentially the same 
as those shown in Fig. 1, that is, q = 1.055 ± 0.002. 

Then, we investigated the possibility that the discrepancies 
observed between WMAP data and exact Gaussian distribu¬ 
tions, as evidenced in Figs. 1 and 2, occur just by chance. 
For this scope, we analysed a set of 10, 000 Monte Carlo re¬ 
alizations of CMB Gaussian maps, and we found, at the 99% 
confidence level (CL), that the distributions of CMB tem¬ 
perature fluctuations measured by WMAP are not properly 
described by exact Gaussian temperature distributions. 

We also investigated, through substantial numerical simu¬ 
lations, the possibility that the instrumental noise could in¬ 
troduce non-Gaussian signals in the CMB maps. According 
to instrument noise produces a random Gaussian sig¬ 

nal, actually Gaussian per observation. Since pixels are ob¬ 
served a different number of times, the effect on the CMB 
map could be a significant non-Gaussian signal in the CMB 
temperature distribution, and clearly this possibility should 
be taken into account in the simulations. For this, to each one 
of the 10, 000 Monte Carlo CMB Gaussian maps we added a 
simulated non-stationary Gaussian radiometer noise, taking 
into account the actual number of observations ( rn ) for each 
pixel in the maps. We consider the rn data from the WMAP- 
W4 map jfi]. We obtained q = 1.005 ± 0.01 for the CMB 
plus non-stationary Gaussian noise simulated maps, which is 
obviously consistent with Gaussianity. A \ 2 estimator test 
shows that the P q (AT) does not fit the simulated data as 
well as it does the WMAP data. The monotonic discrepancy 
behavior observed in Figs, lb, Id, and If, which is well fit¬ 
ted by the nonextensive expression given in eq. 0 due to a 
minimum \ 2 value, was observed in less than 1% of the sim¬ 
ulations performed. Therefore we can rule out, at the 99% 
CL, the simulated radiometer noise as being the explanation 
for the non-Gaussianity observed in the WMAP temperature 
fluctuations. 

In conclusion, we have shown that an exact Gaussian dis¬ 
tribution is excluded, at the 99% CL, to properly represent 
the CMB temperature fluctuations measured by WMAP. Al¬ 
though the value of the g-parameter is close to 1, our analy¬ 
ses indicate that to consider these temperature fluctuations 
as being of Gaussian nature is not rigourously exact and it 
should be considered as a good approximation instead. 
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FIG. 1: Fits to CMB temperature fluctuations measured by WMAP in bands Q, V, and W (after using the KpO mask) 
applying an exact Gaussian distribution (blue curve) and a nonextensive function (red curve). In figures (a), (c), and (e) we 
plotted the normalized number of pixels versus AT 1 , while in figures (b), (d), and (f) we plotted the normalized number of pixels 
versus (AT/ov) 2 , respectively. The \ 2 best-fit for the distributions P q gives q = 1.056, with ctq = 140.58/rK, cry = 163.52/rK, 
and aw = 190.35/rK, respectively. 
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FIG. 2: Results of the analysis of WMAP data using the individual radiometer maps Ql, Q2, VI, V2, Wl, W2, W3, and W4, 
after applying the KpO mask. In all plots the blue curve represents the best-fit Gaussian, according to the estimator, while 
the red curve represents the best nonextensive distribution function. Each P q distributions fits the corresponding data with a 
slightly different value of q, where the mean value of them is q = 1.055 ± 0.002. 
















